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In numerically constructing a spacetime that has an approximate timelike Killing vector, it is 
useful to choose spacetime coordinates adapted to the symmetry, so that the metric and matter 
variables vary only slowly with time in these coordinates. In particular, this is a crucial issue 
in numerically calculating a binary black hole inspiral. An approximate homothetic vector plays 
a role in critical gravitational collapse. We summarize old and new suggestions for finding such 
coordinates from a general point of view. We then test some of these in various toy models with 
spherical symmetry, including critical fluid collapse and critical scalar field collapse. 



Q ■ I. INTRODUCTION 

o; 

^ One of the fundamental problems in numerically constructing spacetimes without symmetries (3D numerical rel- 
• ativity) is the choice of suitable coordinates. The task is not just to find a good coordinate system for a known 
spacetime, but the spacetime is not actually known when the calculation begins, so that one has to construct the 
^sj . coordinate system along with the spacetime. In principle, once one has obtained the spacetime in some coordinate 
^ ' system, one can go back over it and paint on a better coordinate system. In practice, the initial attempt must not be 
, too bad, or the coordinate system will break down altogether before one has evolved very far. 

One begins the numerical construction of a spacetime with Cauchy data (three-metric gah, extrinsic curvature Kab, 
and perhaps suitable matter data) obeying the Hamiltonian and momentum constraints, in hand. One then has to 
■ make a choice of the lapse a and shift based on this information in order to start constructing the spacetime along 
I with the coordinate system. For the purpose of this paper, a coordinate condition is a prescription that maps Cauchy 
data on a Cauchy surface to lapse and shift fields on that Cauchy surface. One then uses the same prescription again 
^ ' on each time slice. We call such coordinate conditions local in time. Typically they are elliptic (in space), and one 
^H, specifies suitable boundary conditions at the outer (large radius) and inner (black hole excision) boundaries of the 
I ' numerical domain. 

In a generalization of local-in-time coordinate conditions, one may use a prescription for the first or second time 
derivatives of the lapse and/or shift, obtaining a heat equation or wave equation-like prescription. Such evolution 
equations for the lapse and shift require much less computational work than solving an elliptic equation at each time 
step. Coordinate conditions for the purpose of our paper are then elliptic, parabolic or wave equations for the lapse 
^ ' and shift in terms of given Cauchy data on a slice. 

What are the characteristics of good coordinate conditions? Here we propose the following desirable criterion: 
If the spacetime is in fact stationary, the coordinate condition should make the metric coefficients explicitly time- 
independent, independently of how the initial slice is embedded in the spacetime. If there is an approximate timelike 
Killing vector, the metric should be approximately time- independent. If there is no Killing vector, the natural 
extension of our proposal is to minimize some measure of the rate of change of the metric with time. 

Consider the inspiral phase of a binary system. There is an approximate timelike Killing vector here, in the sense 
that all metric coefficients can be made to evolve on the inspiral time scale instead of the much shorter orbital time 
scale. Clearly, this defines a corotating coordinate system. If the spacetime has an approximate periodicity (discrete 
isometry connecting timelike related points), as in the case of a highly elliptic inspiral, the metric should become 
approximately periodic in time through our coordinate condition. 

A coordinate system adapted to the presence of a timelike Killing vector, as we have just described it, is far from 
unique. It is easy to see that the family of such coordinate systems is equivalent to the family of smooth spacelike 
slices containing one arbitrary but fixed spacetime point, together with spatial coordinates on the slice: the coordinate 
system on the spacetime is obtained by Lie-dragging the slice and its coordinates along the Killing vector. 

One may know that the final state of an evolution (for example a binary black hole system) is a Kerr solution, but 
this does not mean that the metric coefficients approach the Kerr metric in any of the standard coordinate systems 
even as they become asymptotically time- independent. We require only that our coordinate condition evolves any 
given initial slice along the timelike Killing vector, not that it "straightens out" the slice. A "straightest" slice in this 
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spacetime could be defined as one that minimizes for example J R dV, but we sfiall not try to minimize this at the 
same time. 

A mathematically similar problem to that of finding good coordinates in numerical relativity arises in the theory of 
critical phenomena in gravitational collapse. The entire rich phenomenology of critical collapse can be explained and 
quantitatively predicted using ideas taken from renormalisation group theory and the theory of dynamical systems 
Implicit in this theoretical framework is the concept of a coordinate prescription (local in time) such that the 
metric is "time" -independent if the Cauchy data generate a spacetime that is homothetic, or continuously self-similar 
(CSS). We have set "time" in quotation marks here, as d/dt is not everywhere timelike in a homothetic spacetime. 
In this context we are only interested in a prescription that is local in time (in the sense of not containing time 
derivatives of the lapse and shift). Such a prescription turns general relativity into a dynamical system, with CSS 
spacetimes corresponding to fixed points. If the spacetime is discretely self-similar (DSS), we require our prescription 
to make the metric periodic in "time" , with a DSS spacetime corresponding to a limit cycle. (The use of maximal 
slicing in critical gravitational collapse Q does part of the job of turning general relativity into the appropriate sort 
of dynamical system. However, one is still left with the question of what spatial coordinates to use). 

Again, such a coordinate system is far from unique. In the neighborhood of any data set generating a CSS or 
DSS spacetime one can define such a coordinate prescription perturbatively. This perturbative definition has been 
enough in order to calculate critical exponents to several significant digits. A coordinate condition applicable to all of 
superspace, with the property that the fixed points of the resulting dynamical system are precisely CSS spacetimes 
and the limit cycles the DSS spacetimes, would define a renormalisation group for general relativity. Studying its 
properties would help us in understanding the general dynamics of general relativity and issues relating to cosmic 
censorship. 

In this paper, we study the two related cases in a common framework. In section II, we discuss a number of 
candidate coordinate conditions in the absence of symmetries (other than stationarity or homothety). Our basic 
strategy is to write down the condition that the metric is (conformally) time-independent in a 3-1-1 split, and then 
choose four equations out of that set of equations to solve for the lapse and shift or their time derivatives. This will 
guarantee that our criterion is obeyed. If a given prescription obeying this criterion is otherwise a good coordinate 
condition must be determined empirically. In section III, we specialize the general proposals to spherical symmetry, 
and we discuss additional proposals that are specific to spherical symmetry. In section IV, we test some of these 
proposals by evolving a slice in a given spherically symmetric spacetime. (This procedure is simpler and more stable 
than evolving the Cauchy data themselves.) We use several spacetimes to test the stationary case, the CSS case, and 
the DSS case. 



II. COORDINATE CHOICES 



A. 3+1 split of the Killing and homothetic equation 



We begin with a homothetic Killing vector X"'. That is 

jO-xgab = Va^fc + VfoXa = 2(7gab, (1) 

where cr is a constant. For a — 0, A"" is a Killing vector. For a ^ a. constant we can, without loss of generality, 
choose (7 = — 1 by rescaling X". In the following a only takes the values and —1. We perform the usual 3-1-1 split 
of the Einstein equation with X" as the evolution vector field. That is, spacetime is foliated by surfaces 'S{t) with 
AT" ~ {d/dty\ Each surface has a unit normal vector n", intrinsic metric hab — gab + nani, with associated derivative 
operator Da and extrinsic curvature Kab — —{^/'2.)Cnhab- (There is no standard convention for the sign of Kab- Here 
we use the sign prevalent in numerical relativity.) We decompose X"" as art" -I- with (3°'na = 0. We want to know 
how the quantities {a, (3°^ , hab, Kab) evolve with time. We will get one answer from the homothetic Killing equation 
and another answer from the usual ADM equations. The condition that these two answers agree then gives conditions 
that these quantities must satisfy in a spacetime with a homothetic Killing vector. . 

We begin by deriving evolution equations from the homothetic Killing equation. From the decomposition of AT", it 
follows that Ua = —aS/at and therefore that 

Cxria = a^^{Cxa)na. (2) 
Since Cx^a is a scalar times Ua, it follows from equation (|l]) that 

Cxria = ana, (3) 
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and therefore that 

CxOL — aa. (4) 

Then using = ^x^" = we find 

CxP" = 0. (5) 

It follows from equations (0) and (^) that 

Cxhab = 2ahab- (6) 

We then find 

^xKab — —-CxCnhab — — {^n^xhab + ^[X,n]hab) , (7) 

and therefore 

CxKab = aKab. (8) 



B. Solving for the lapse and shift locally in time 

We begin by looking at those components of the homothetic Killing equation that contain time derivatives of the 
lapse and shift. With the shorthand CxT = T for any tensor T, these are 

a = aa, = 0. (9) 

Given a good initial guess, we could consider these as evolution equations for the lapse and shift, but this is unlikely to 
make for a stable scheme when the symmetry is only approximate and/or the initial guess is not perfect. Furthermore, 
it leaves us with the problem of finding a good initial value for the lapse and shift. Conversely, if we have found a lapse 
and shift on the initial time slice by solving some elliptic equations, the only reason not to use the same prescription 
on each time slice would be numerical effort. [If numerical effort is a concern, one could in effect evolve the lapse and 
shift by (^) for a few time steps, before using the expensive prescription again.] 

We now concentrate on the other components of the homothetic Killing equations. The ADM evolution equations 
are 



Kab = P'DcKab + KacDbP' + KbcDafi" " DaDbU + a I ^'^^Rab + KKab - 2KacKh' - Tab - \{p - T)hab ] , (11) 



DaPb + DbPa - 2aKab, (10) 

1 

2< 

for any vector field X°-. Here p = Gabn'^n'', Tab = ha'^hb'^Gab, and Ja = —ha'^n^Gcb- Note that in the ADM evolution 
equations, a and /J" can be freely specified, so no new information is obtained by combining the ADM equations with 
the equations (||). But if we substitute the ADM equations into the components 

hab = 'i.ahab (12) 

kab = CjKab. (13) 

of the homothetic Killing equation, we obtain two symmetric tensor equations which link the lapse and shift, and 
their spatial derivatives, to the Cauchy data, and their spatial derivatives, all within a single Cauchy surface. 

Our task can now be formulated as follows. Out of these 6+6 equations we want to obtain 3+1 equations that can 
be solved for the lapse and shift, for given generic Cauchy data. The remaining 8 equations remain as consistency 
conditions which are obeyed if and only if the Cauchy data actually evolve into a spacetime with a (conformal) Killing 
vector. In other words, given Cauchy data on a slice, we want to have a prescription to calculate a lapse and shift 
such that one obtains a reasonable coordinate system on the spacetime evolved from "reasonable" generic data, and 
coordinates adapted to the symmetry if there is a timelike Killing vector or homothetic vector. 
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Clearly, there is an infinite range of possibilities. In a situation without symmetries, our task is to make one 3- vector 
(the equation for the shift) and one scalar (the equation for the lapse) out of two symmetric 3-tensors. It is important 
to note that the equations for the lapse and shift will generally be coupled, so that the same equation for, say, the 
shift can have a different character when coupled to a different equation for the lapse. Some or all of the resulting 
equations can be purely algebraic. Alternatively, they can involve first or higher spatial derivatives of the lapse and 
shift. An elliptic equation or equations would be particularly appealing. Any prescription obtained for the general 
case can be reduced to spherical symmetry, and we shall do this for testing. Additional prescriptions that explicitly 
use the spherical symmetry do not generalize back to the case without symmetries. 

One can obtain a 3-vector of equations from a 3-tensor by taking a divergence, for example 



D'' [hab - 2ahab) = D^hab = 0. (14) 

This equation can be derived from varying the action 

J d^xVhh''^h'"'ih,b - 2ahab){hcd - 2ah,d) (15) 
with respect to the shift. As a drops out of (^, it could also have been derived from varying the action 

(f'x\/hh'"'h^'^habhcd- (16) 



Condition ( p^ was first suggested as a shift condition by Smarr and York Q under the name minimal strain shift. 
As it is derived from a variational principle, one may hope that it gives good spatial coordinates also in the absence 
of an exact symmetry. Written out explicitly, minimal strain shift is 

DbD^P"" + DbD^p^ - 2Db{aK''^) ^ 0. (17) 

Taking the divergence of (|l^), we obtain the 3-vector equation 

D'' (kab - CjKab) = 0. (18) 

At first sight this is less appealing than (|lj) because it contains third spatial derivatives of the lapse. 
In order to obtain a scalar equation for the lapse, we can take the trace of (p^; 

(hab - 2ahab) = (19) 

is equivalent to 

Vh-3(7Vh = 0, (20) 

that is, it determines the scaling of the 3-volume element, {h is the determinant of hab-) Written out in terms of a 
and this equation is 

DaP" - aK - 3a = 0. (21) 

This does not seem to be in general a good equation for the lapse (since it is ill defined where K vanishes), nor (by 
virtue of being a scalar) for the shift. 

The equation resulting from instead contracting (|lj) with Kab, 



(hab - 2(jhab) - (22) 



K 

can be derived from varying the functional ( p^ ) with respect to the lapse. We shall refer to Eq. ( p^ ) on its own as 
the conformal strain lapse for cr = —1, or the minimal strain lapse for cr = 0. It is an algebraic equation for the lapse, 

K'^^DaPb - (tK 

which is linear in the shift. The fact that the lapse may be zero, negative, or even blow up (where Kab = 0), is 
worrisome. Clearly this prescription cannot be used on all initial data. For the Killing case, c = 0, Brady, Creighton 
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and Thorne (BCT) |^ have recently proposed using the minimal strain lapse and minimal strain shift together, as 
both are derived from the same variational principle (^^. If one substitutes the formal minimal strain lapse (|2^) into 
the minimal strain shift equation (p7|), it is not clear if the resulting differential equation for the shift alone is elliptic. 
We shall refer to the combination of minimal strain lapse and minimal strain shift, in the Killing case ct = 0, as the 
BCT gauge, and for any a as the generalized BCT gauge. 
The equation 



aKab]=Q (24) 



is an elliptic equation for the lapse, and seems not to have been considered before. We shall call it the trace-Kdot 
lapse. Adding to this the minimal strain lapse equation, however, we obtain 

- h'^'' (kab - <jKab) - K"'' (hab - 2ahab) =k + aK. (25) 



For (7 = 0, this is known as constant (extrinsic) curvature slicing. If if = initially, and therefore at all times, we have 
maximal slicing. For cr = — 1 we could speak of exponential curvature slicing. Written out in terms of the Cauchy 
data, the equation is 



P^DaK + aK ^0. (26) 



Note that we do not assume that K is constant within each slice. We shall refer to the combination of equations (|ri_ 
and (|2^) as the generalized Smarr-York (SY) gauge. For an initial K — slice, the generalized SY gauge reduces to 
maximal slicing with minimal strain. 

The differential equations for lapse and shift must be supplemented with boundary conditions in order to yield 
a solution. For an approximate timelike Killing field in an asymptotically flat spacetime, the condition a — > 1 and 
Z?" ^ at infinity should be reasonable. For the case of binary black hole inspiral, BCT have proposed a boundary 
condition to use with their equation. It is less clear what would be a reasonable boundary condition to use for the 
case of critical gravitational collapse. 



III. SPHERICAL SYMMETRY 



A. Reduction to spherical symmetry 



In spherical symmetry, the metric has the form 

ds^ = (-a^ + (3'^(3r)dt^ + 2/3rdrdt + grrdr"^ + gee{dO'^ + sm^Odcj)^) 



(27) 



where a,l3^ ,grr and ggg are functions of r and t, and where Pr = grrP^ ■ Equations (|l^ ) and ( |l3|) now have only two 
independent components each. The rr and 96 components of equation (03) are respectively 



1 d-drr 

2' dr 



dr 



- crgee = 0. 



l^' dr 



The rr and 96 components of equation (nSf) are respectively 



1 dgrr da ^ dK, 
P 



dr"^ 2grr dr dr 



dr 



2Kr 



'^rr , Qrr'^OO P9rr 



gee grr J 2 



1 dgee da ^ dK, 



2grr dr dr 
The constraint equations 



dr 



(3) 



KrrKee f ]jr 

grr 2 y g^ 



dr 

aKrr — 0, 
<jKee = 0. 



(28) 
(29) 



(30) 
(31) 
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+ 2 _ j^^^^j^ab _ 2p = 0, DaK\ - D^K - Jfc = 0, 

become in spherical symmetry 

P = -7y + 

^9rr 900 
_ dKf)B 

—geeJr 



dr 





f2Krr 


^Kee\ 


900 


V 9rr 


900 J 




^ Kee\ 


dgee 


\ 9rr 


960 ) 


dr ■ 



(32) 

(33) 
(34) 



B. Generic coordinate conditions applied to spherical symmetry 



In spherical symmetry, the conformal strain lapse, equation (23) is given by 



2 Kee 

gee 



Srr / \ gee 



The exponential curvature lapse, equation (p 



is given by 
dr^ 



+agrr 



+9r 



Krr 
9rr 

d /Kr- 



2(^] 
\ 900 ) 

Ke 



\2g rr 

2 



1 dgee 

dr gee dr 



da 
dr 



grr 



Tee 



g00 



dr V grr gee 
The minimal strain shift, equation (0^), is given by 



Kr 



d'^p'- 



+- 



1 d^grr 
Qrr dr^ 



1 dgee 
gee dr 



dr^ 
1 dgr, 



dr 



,Kee 

' g00 



1 dgr 



= 0. 



+ 



1 



2grr dr ge0 dr 



gr 



dr 



1 dgrr 

Qrr dr 



1 dgee 
100 dr 



\- 


( aKrr\ 


1 a dgee 


( K 


Kee\ 


dr 


K grr ) 


gee dr 


\ 9rr 


900 ) _ 



= 0. 



(35) 



(36) 



Thus, in spherical symmetry the generalized SY gau ge i s the combination of equations ( p6[ ) and (^ 
generalized BCT gauge is the combination of equations (p^) and (^7|). 



(37) 
while the 



C. Coordinate conditions tailored to spherical symmetry 

In spherical symmetry, one can also find simpler equations for the lapse and the shift. We note that equation ( [29| ) 
contains no derivatives of the shift, while equation (|3l|) contains only one derivative of the lapse. We can use equation 
( p9| ) to express the shift as 

.(^900 + aKe0 



(3'' = 2- 



(38) 



dge0/dr 

For a — Q this is called the area freezing shift, as it keeps the area of a surface for constant r constant. Similarly, we 
call ( |3l| ) the area lapse. Substituting ( pq ) into (pl|), we obtain a first-order differential equation for a alone. We can 
use the momentum and Hamiltonian constraints (p4, |33[) to simplify it to the form 



da 

dr 



dgee/dr 



2grrKeeJr 



+ Tr 



dgee/dr J \ dgee/dr 

We can integrate this first-order ordinary differential equation from the center of spherical symmetry outwards. 
However, the equation requires as a boundary condition the value of the lapse at one value of r. We refer to the 
combination of equations ( p8| ) and ( |39| ) as the area gauge. 



2a K, 



grrgeeJr 



(39) 
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D. Boundary conditions for the lapse and shift equations in spherical symmetry 



In three space-dimensions without spherical symmetry, some of the equations for the lapse and shift we have dis- 
cussed are coupled elliptic equations. These require boundary values for the lapse and shift on a surface (topologically 
a 2-sphere) at large radius. One would expect the equivalent problem in spherical symmetry to be a second-order 
ODE boundary value problem, with one boundary condition each at r = and large r, for each of a and (3"^. One 
boundary condition at r = is immediately obtained from symmetry, namely that a must have an expansion in even 
powers of r and an expansion in odd powers of r. For the second boundary condition, it is as natural to specify a 
and lim(/3'"/r) at r = as it would be to specify a and (3"^ at some large value of r. Doing that, we have the great 
simplification of integrating an ODE system from r — outwards, instead of solving a boundary value problem. 

The equations for a and are homogeneously linear in a, and a. The boundary values we need can be thought 
of as an overall factor in both a and P^ /r, and their ratio. The ratio is uniquely determined at the center of spherical 
symmetry. In general we have 



aK + DaP" 



9(lnVg) 

dt 



3ct, 



(40) 



where the first equality follows from the definition of K, while the second equality follows from (^^. At the center of 
spherical symmetry, this yields 



P'' 



a + iisTa ) r + 0[r^ 



(41) 



The same result is also obtained by evaluating either (|2l|), ( p3[ ) or ( |38| ) at r = 0. 

It remains to determine the overall factor in both a and /?''. From (^) we have that a = ao = const, at the origin 
in the Killing case. The value of ao is a matter of convention that only affects the norm of the Killing vector field, 
and we can set it to any value that is constant in i. 

In the self-similar case, we formally have a = ao exp at at the center, but this result is not very useful, as only 
one choice of olq gives rises to a compatible coordinate system. To see this, label the central geodesic by proper time 
T with the origin chosen so that T = is the singular point of the self-similar spacetime. In the self-similar time 
coordinate t this point is t = oo. Let T — Tq where the initial slice t = intersects the central geodesic. Then it is 
easy to see that the proper relation T — Toexpat is obtained only if we choose aa — (tTq. In other words, we must 
guess correctly how far the slice is away from the singular point. 

A prescription for a at the center can be obtained by noting that we have CxTab — and therefore CxP = —'lap. 
In the 3-f 1 split this is 



We combine this with the equation for n'^\I aP obtained from stress-energy conservation, 

= -UaVbG"^ = -Kp -f n^VaP + - T'^'^Kab + ^J^Da In a. 



(42) 



(43) 



At the center of spherical symmetry, DaCi and P"' both vanish, and we obtain an algebraic equation for a alone. 
Furthermore at the center we have Kab = {K/3)gab and a similar expression for Tab, which simplifies our final result 
to 



2apgr 



dr 



rr \ P+ — 



+ 0(r2) 



(44) 



This equation gives a (nonzero) value for a at the origin provided that neither the numerator nor the denominator 
vanishes there. 



IV. TESTS IN SPHERICAL SYMMETRY 



A. Numerical method 



We now discuss empirical tests of some of the coordinate conditions we have discussed for some given spacetimes, 
restricting ourselves to spherically symmetric examples. Rather than to evolve Cauchy data via the ADM equations, 
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we evolve an embedded slice in the spacetime. All indices A, B, C, i,j, k in this section denote coordinate components 
of tensors, rather than abstract tensor indices. Let a spacetime metric be given in closed form in coordinates x^, 
giving rise to a metric gAB- A three-dimensional spacelike slice embedded in the spacetime, together with a choice 
of three coordinates on the slice, is parameterized by four functions of three coordinates, = F^{x'^). The 
four-metric qab induces a three metric gij and extrinsic curvature Kij , and there is a unit timelike normal to the 
slice at each point. 

In our code, we choose a spacetime, and an initial slice, and calculate and Kij. The equations for doing this are 
derived in full in If the spacetime contains matter, the components p, Ji and Sij of the stress-energy tensor are 
also computed from Tab- These are handed to a coordinate condition solver that returns a and The slice is then 
evolved via 

We finite-difference the equations using a type of iterative Crank-Nicholson (ICN) algorithm. dF^ /dt is computed at 
time tn, and the are evolved to i„+i with a forward in time, centered in space (predictor) step. dF^ jdt is then 
computed at time and the average of dF^ /dt at t„ and tn+i is then used for a corrector step. The corrector 

step is iterated to convergence. To maintain stability, we use numerical viscosity. 

In spherical symmetry, there are only two functions X''\r) and X^{r), and two components each of the induced 
metric and extrinsic curvature, grrir), geeir), K^rir), and Kgg{r), and stress-energy components p{r), Jr(r) and 
Trrir). There is only one shift component, f3^(r), and the lapse, a{r). The coordinate condition solver and slice 
evolution code are independent and communicate only through the fields intrinsic to the slice, so that the coordinate 
condition solver never sees the F^ directly. The code assumes that surfaces of constant t are spacelike, and complains 
if, through a bad choice of coordinate, they are not, but no assumptions are made on the coordinates X° and X-^. 
They could, for example, be null coordinates. 

Integrating the embedding of a slice, rather than Cauchy data, has two advantages. The first one is stability: 
empirically, the ICN algorithm, together with a linear extrapolation outer boundary condition, is very stable for given 
a and . In particular, if we want to excise the center of spherical symmetry from the numerical domain in a black 
hole spacetime, an extrapolation inner boundary works well in the slice evolution code, while black hole excision in a 
genuine spacetime evolution code is a nontrivial problem even in spherical symmetry. The second advantage of the 
slice evolution is that we can use a given spacetime that is not a solution of the Einstein equations for any reasonable 
matter, but that has some properties worth investigating. In such a spacetime, we create a fictitious stress-energy 
tensor simply defined by the Einstein equations but not derived from any reasonable matter source. We think that 
such "fake" spacetimes can be valuable tests of qualitative aspects of coordinate conditions. 



B. Stationary spacetimes 

We have used a variety of spacetimes to test three different gauges: area gauge, generalized SY gauge and generalized 
BCT gauge. The results are as follows. 

As model stationary spacetimes with a regular center we have examined Minkowski space, de Sitter space and 
the Schwarzschild solution for a constant density star, parameterized by -RstarZ-^star < 9/4. All three gauge choices 
work on the Minkowski and de Sitter spacetimes. The area and generalized SY gauges work on the constant density 
star. The generalized BCT gauge fails numerically at the boundary of the star, where dgrr/dr is discontinuous. As 
a boundary condition we choose a = 1 at the center. As the equations are linear in both a and (3^, this choice is 
equivalent to any other choice that is constant in time. As a consistency check of our programming, we have verified 
that the Hamiltonian and momentum constraints are obeyed up to finite differencing error. As the initial slice, we take 
X^ = r and X^{r) a Gaussian, with X° and X^ the standard time and radial coordinates. We find that the induced 
metric, extrinsic curvature, stress-energy components, lapse and shift are all time-independent up to numerical error. 

To test coordinate conditions in a black hole spacetime with the black hole excised from the numerical domain, 
we have implemented the Schwarzschild spacetime in coordinates that are regular at the future event horizon (the 
Painleve-Giillstrand coordinates). Here the area gauge and generalized SY gauge work. The generalized BCT gauge 
has numerical problems at the excision boundary. As a boundary condition we have chosen a = 1 at the excision 
boundary, which is chosen at constant coordinate r. As the shift is supposed to keep our generic radial coordinate 
r at constant curvature radius ^/gee, this boundary condition is consistent. As the initial slice, we took X^ — r and 
X'^{r) a Gaussian, with X^ the area radius and AT" the Painleve-Giillstrand time (which is defined by surfaces of 
constant A° being flat and spacelike). Again we find that all variables are constant in time, even though they are all 
nontrivial. 
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C. Self-similar spacetimes 



We have also investigated four self-similar spacetimes. Our coordinate prescriptions are designed to bring the metric 
of a CSS spacetime into the form 

5^,. = e-2*g^^ (46) 

where the components of are independent of t. The metric of a DSS spacetime can be written in the same form, 
with g^t^ periodic in t. It is important to note that while all our prescriptions must work on CSS spacetimes, barring 
numerical error and numerical instability, it is not guaranteed that they work, in the sense of bringing the metric into 
this form, on DSS spacetimes. The reason is that given a CSS spacetime and an initial slice, the coordinate system 
of this form is unique, while it is not for a DSS spacetime. In order to see how successful a gauge is in making the 
conformal spacetime metric independent oit or periodice in t, it is useful to rescale fields intrinsic to the slice (gij, Kij 
etc.) with appropriate powers of expi so that the rescaled fields are independent of time if and only if the spacetime 
is CSS and our algorithm finds the conformal Killing vector. 

The simplest CSS spacetime, Minkowski spacetime, is not suitable for our investigation, as the homothetic vector 
is not unique. Instead, we have considered a conformally flat spacetime of the form 

We have given this spacetime fictitious matter defined by the Einstein equations. It has a unique homothetic vector 
field if by assumption we restrict that vector field to the tr plane. Here, we have tested all three gauges. All work 
well when we use the the analytic boundary condition a = e~' at the center. 

A less trivial example is the Roberts [|| solution, a family of CSS massless scalar field solutions, which in double 
null coordinates is given by 

ds^ = ~du dv+^ [(1 - p^)u^ - 2uv + u^] dn^ (48) 

cp{u,v)^ln\—^ , (49) 

(1 + p)v — u 

with p a parameter. This spacetime has curvature singularities at gee = 0. We have investigated the case p = 0.5, 
where these singularities are naked, and lie in the future and past lightcones. The slice we evolve lies in the spacelike 
sector. The only boundary condition for a at the inner boundary that worked here is a = ao exp at, with ao taken 
from the exact solution. Here the area gauge and generalized SY gauge work. The generalized BCT gauge has 
numerical problems at the excision boundary. 

Our other two examples are of more physical interest. They are the CSS critical solution for the collapse of a perfect 
fiuid with the equation of state p = (l/3)p and the DSS critical solution for the collapse of a massless scalar field 
. These solutions have been constructed numerically as solutions of a nonlinear ODE eigenvalue problem |^ , and 
a nonlinear PDE eigenvalue problem , respectively. We interpolate them numerically to the required values of 
and X^. For both spacetimes, we have tested the area gauge and the generalized SY gauge. Again, we use a Gaussian 
initial slice. 

In the CSS fluid spacetime, we have been able to use the boundary condition (^). The boundary condition works, 
and all variables are constant in time. 

The DSS scalar field spacetime has a regular center, but the scalar field is a periodic function of r at the center, so 
that the density p vanishes twice per period. This means that we cannot apply the boundary condition (^4|). Instead, 
we have applied the exact boundary condition a — ao expert, with ao taken from the exact solution. This test is 
nontrivial, as it is not clear at all that a coordinate condition designed to catch a continuous symmetry can catch a 
discrete one. But it works: all variables are periodic in time for the scalar field DSS spacetime for at least two periods, 
in both the area and generalized SY gauge. 



V. CONCLUSIONS 



We have examined the problem of finding good coordinates in numerical relativity, and the problem of defining a 
renormalisation group for general relativity in critical collapse, in a single framework. We have listed old and new 
prescriptions that obey our criterion for a good coordinate condition in either application: that the presence of an 
approximate timelike Killing field, or an approximate homothetic vector field, in some region of the spacetime (for 
example at late times) be revealed by metric coefficients that vary only slowly with time. 
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We have tested several prescriptions numerically in spherically symmetric spacetimes that are exactly stationary 
or homothetic. For these exact symmetries, it is clear analytically that our prescriptions will work, but we have used 
these cases both to test our code, and to see how well different schemes can be implemented numerically. We find 
that a scheme specific to spherical symmetry (area freezing shift and radial lapse) is the most stable one, generalized 
Smarr-York (SY) gauge (minimal strain shift and exponential K lapse) comes second, while we have been unable to 
implement the generalized Brady-Crighton-Thornc (BCT) gauge (minimal strain shift and conformal strain lapse) 
stably. What this means is that the combination of our numerical methods of solving for the BCT gauge and of 
evolving a slice is numerically unstable. This indicates that coupling the BCT gauge to a genuine spacetime evolution 
code may also be numerically unstable, but it is of course possible that a different numerical implementation of the 
gauge and/or the evolution equations is stable. Other difficulties of the BCT gauge are evident analytically: Not all 
initial data sets (in our case, initial slices through a given spacetime) admit a solution to BCT gauge, and BCT gauge 
may still break down (when KabK"''' = anywhere) at a later time. In spherical symmetry, the slice must nowhere 
be close to a polar slice for BCT gauge to work. 

We have entered new terrain examining one spacetime (the Choptuik scalar field critical solution) with a discrete, 
rather than continuous, symmetry. Here we found numerically that two prescriptions made for a continuous symmetry 
(area gauge and generalized SY gauge) actually follow the discrete symmetry, in the sense that the metric becomes 
periodic in time. (Our third gauge, generalized BCT gauge, was unstable numerically in our implementation, but 
may yet work with better numerical methods.) This is an important step in the construction of a renormalisation 
group for critical collapse, as it suggests that a single coordinate prescription may realize the renormalisation group 
both for continuously and discretely self-similar (CSS and DSS) critical solutions. It is also encouraging for the use 
of co-rotating coordinates in non-circular binary black hole inspiral. 

One crucial ingredient is still lacking in a working prescription for a renormalisation group, namely the boundary 
conditions used when solving elliptic equations for the lapse and shift. Details aside, these boundary conditions must 
contain a good guess for how far away in proper time the slice is from the accumulation point of the (approximately) 
self-similar spacetime. Although there are ideas [see equation (^)], we have not been able to solve this problem, and 
instead had to supply this information by hand. 

Future work should include testing beyond spherical symmetry and in particular with angular momentum, and on 
spacetimes with approximate symmetries, as well as developing ideas on what boundary conditions should be imposed 
on elliptic coordinate conditions. 
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